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A computer  program  to  compute  estimated  orbital  element 
sets  based  on  positional  observations  from  a single  observing 
site  has  been  developed  and  tested.  It  provides  a reasonable 


mix  of  operator  interaction  and  automation,  making  it  adaptable 
to  the  level  of  proficiency  of  the  user,  r Included  here  are'a 
program  overview,  concept  of  operations,  specific  operating 
instruction,  examples  of  program  operation  and  examples  of 
program  output. 
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OVERVIEW 


I . 

The  initial  orbit  estimation  module  (IOMOD)  is  an  interac- 
tive task  which  allows  the  user  to  generate  classical  orbital 
element  sets  from  satellite  positional  measurements.  The  input 
data  consists  of  angles  in  either  the  right  ascension/declination 
or  azimuth/elevation  reference  frame,  and  optional  range.  These 
data  may  be  read  from  the  system  observation  file  or  from  the 
card  reader  or  terminal  keyboard. 

The  program  contains  an  algorithm  for  range  estimation  from 
angles-only  data  and  two  orbit  estimation  algorithms.  Range 
estimation  is  via  the  classic  Gauss-Gibbs  technique.  The  primary 
orbit  estimation  method  is  a modified  true  anomaly  iteration 
technique.  This  method  historically  produces  usable  orbit  esti- 
mates for  the  classes  of  satellites  of  interest  to  the  ETS. 

Being  iterative,  it  can  occasionally  fail  to  converge.  For  these 
times,  and  for  an  independent  consistency  check,  the  three  obser- 
vation Herrick-Gibbs  technique  is  supplied. 

Default  values  for  the  sensor  coordinates  and  default  pro- 
cessing methods  are  available.  The  inexperienced  user  can  rely 
on  the  program  to  produce  results.  The  more  experienced  user 
will  find  sufficient  flexibility  to  put  his  expertise  to  use. 
Among  the  features  available  are:  buffering  of  up  to  twenty 
observation  sets;  editing  of  the  observations  (ADD/DELete/ 
REPlace/change) ; observation  selection  for  use  in  orbit  estima- 


1 


tion  method;  satellite  tagging  and  orbital  element  set  file 
up  capability. 


V 


II.  OPERATION 

IOMOD  is  available  to  both  the  batch  user  and  the  real-time 
operator.  The  batch  version  is  cataloged  as  IOD  (Initial  Orbit 
Determination) , while  the  real-time  version  is  cataloged  as  RIOD. 
IOD  is  activated  from  the  computer  operator's  console: 

/IOD/ ACT, , TM 

RIOD  is  activated  from  the  real-time  console: 

/RIOD 

Both  versions  are  identical  except  that  RIOD  is  aware  of  its 
real-time  roll.  More  details  are  provided  in  the  Programmer's 
Notes . 

Once  activated,  the  program  erases  the  CRT  screen,  writes 
its  name,  version  ID,  the  current  GMT  day/date/time  and  the 
default  sensor  coordinates  (Stallion  Range  Center,  A dome),  a 
prompting  question  mark  and  awaits  user  input.  Typical  program 
operation  involves  reading  observational  data  (OBS  command)  and 
selecting  observations  from  among  those  read  for  orbit  estimation 
(SELect  command) . An  automatic  selection  algorithm  is  available, 
or  specific  observations  may  be  selected  by  line  number.  Each 
SELection  produces  a numbered  CASE  listing.  Summary  output  on 
the  CRT  lists  the  case  number,  observational  data  selected,  iter- 
ative convergence  parameters  and  the  orbital  element  set  produced 
(Figures  1-3  are  typical  terminal  output) . This  and  some  addi- 
tional information  is  sent  to  a printer  file.  One  can  execute 
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INITIAL  ORBIT  ESTIMATION  4.0  DAY  317,  13  Nov  78  12:  4 

SITE,  LAT.  M. , LONG.  E.,  HEIGHT  METERS 
STAL  33.8177  253.3406  1502.4190 

? 

COMMANDS 

+ 

SAT  SIT  OBS  ADD  DEL  REP  R. A DEC  RAN  LIS  SEL  XEQ  YEA  MEE  FIL  XEQ 

? 

YEAR  77 
1977 

? 


83564 

242  19 

242, 

19, 

34 

242, 

19, 

35 

242, 

19, 

36 

242, 

19, 

37 

242, 

19, 

38 

242, 

19, 

47 

242, 

19, 

57 

242, 

20, 

3 

242, 

20, 

17 

242, 

20, 

34 

242, 

20, 

42 

USE 

# 

DAY 

HR 

MIN 

SEC 

RA (DEG) 

DEC (DEG) 

RANGE (KM  OR  ER) 

0 

1 

242 

19 

34 

43.020 

165.84583 

-0.77400000 

0.00000000 

0 

2 

242 

19 

35 

43.020 

166.09583 

-0.77666667 

0.00000000 

0 

3 

242 

19 

36 

37.021 

166.32500 

-0.77833333 

0.00000000 

0 

4 

242 

19 

37 

37.021 

166.57500 

-0.78100000 

0.00000000 

0 

5 

242 

19 

38 

37.022 

166.81167 

-0.78333333 

0.00000000 

0 

6 

242 

19 

47 

13.026 

168.98083 

-0.80833333 

0.00000000 

0 

7 

242 

19 

57 

7.031 

171.47250 

-0.84666667 

0.00000000 

0 

8 

242 

20 

3 

7.035 

172.99833 

-0.87500000 

0.00000000 

0 

9 

242 

20 

17 

41.043 

176.66667 

-0.95433333 

0.00000000 

0 

10 

242 

20 

34 

51.007 

180.98750 

-1.0716667 

0.00000000 

0 

11 

242 

20 

42 

36.012 

182.93750 

-1.1326667 

0.00000000 

? 

SELECT 


83564. 

CASE  1 

STAL 

USE 

# 

DAY 

HR 

MIN 

SEC 

RA  (DEG) 

DEC (DEG) 

RANGE (KM  OR  ER) 

1 

1 

242 

19 

34 

43.019 

165.84583 

-0.77400000 

0.00000000 

1 

8 

242 

20 

3 

7.034 

172.99833 

-0.87500000 

0.00000000 

1 

10 

242 

20 

34 

51.006 

180.98750 

-1.0716667 

0.00000000 

“♦ESTIMATED  MEAN  MOTION  FOR  A CIRCULAR  SATELLITE*** 

83564.0  1.00187268  REV/DAY 

MEAN  SQUARE  RANGE  ERROR  0.273459E-12 

IONU  FAILED 

1 83564U  IOHG3  77242.83549808  0.00000000  0 

2 83564  4.1478  69.5840  9951021  104.7635  358.6330  0.994176614 

? 


Fig.l.  IOD  — console  session,  part  1. 
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SEL 

83564. 


CASE 


STAL 


USE  # DAY 

1 2 242 

1 9 242 

1 11  242 


HR  MIN  SEC 

19  35  43.019 

20  17  41.042 

20  42  36.011 


RA  (DEG) 
166.09583 
176.66667 
182.93750 


DEC (DEG) 
-0.77666667 
-0.95433333 
-1.1326667 


‘“ESTIMATED  MEAN  MOTION  FOR  A CIRCULAR  SATELLITE*** 

83564.0  1.00187268  REV/DAY 


RANGE (KM  OR  ER) 
0.00000000 
0.00000000 
0.00000000 


MEAN  SQUARE  RANGE  ERROR  0. 471657 E- 12 


CONVERGENCE  AT  -0 . 745058E-09 

1 83564U  IONU  77242.86291780  0.00000000  0 

2 83564  4.1542  69.5559  0029964  84.8179  28.3421  0.99930581 

•? 

SEL 


83564 

CASE  3 

STAL 

USE 

# 

DAY 

HR 

MIN 

SEC 

RA(DEG) 

DEC (DEG) 

RANGE (KM  OR  ER) 

1 

3 

242 

19 

36 

37.021 

166.32500 

-0.77833333 

0.00000000 

1 

7 

242 

19 

57 

7.031 

171.47250 

-0.84666667 

0.00000000 

2 

10 

242 

20 

34 

51.006 

180.98750 

-1.0716667 

0.00000000 

** ‘ESTIMATED  MEAN  MOTION  FOR  A CIRCULAR  SATELLITE*** 

83564.0  1.00187268  REV/DAY 

MEAN  SQUARE  RANGE  ERROR  0. 442698 E- 12 


CONVERGENCE  AT  0.909 49 5E- 12 

1 83564U  IONU  77242.85753480  0.00000000  0 

2 83564  4.1480  70.6847  0167755  214.2071  257.9158  1.01615356 

? 

SEL  1,7,9 

83564.  CASE  4 STAL 


USE 

* 

DAY 

HR 

MIN 

SEC 

RA (DEG) 

DEC (DEG) 

RANGE (KM  OR  ER) 

2 

1 

242 

19 

34 

43.019 

165.84583 

-0.77400000 

0.00000000 

2 

7 

242 

19 

57 

7.031 

171.47250 

-0.84666667 

0.00000000 

2 

9 

242 

20 

17 

41.042 

176.66667 

-0.95433333 

0.00000000 

‘“ESTIMATED  MEAN  MOTION  FOR  A CIRCULAR  SATELLITE*** 

83564.0  1.00187268  REV/DAY 

MEAN  SQUARE  RANGE  ERROR  0.301849E-12 


CONVERGENCE  AT  -0 . 102318E- 10 


1 83564U 

2 83564 
? 

XEQ  HG 

1 83564U 

2 83564 


IONU  77242.84561391  0.00000000 
4.1560  71.4156  0294026  219.2969  249.0194 


IOHG3  77242.33133137  0.00000000 
4.1560  71.4156  0294043  219.3011  243.7058 


1.03264931 


1.03265353 


0 


0 


Fig.  2 


IOD  — console  session,  part  2 


LIST 


USE 

# 

DAY 

HR 

MIN 

SEC 

RA (DEG) 

& C (DEG) 

RANGE (KM  OR  ER) 

2 

1 

242 

19 

34 

43.020 

165.84583 

-0.77400000 

0.00000000 

1 

2 

242 

19 

35 

43.020 

166.09583 

-0.77666667 

0.00000000 

1 

3 

242 

19 

36 

37.021 

166.32500 

-0.77833333 

0.00000000 

0 

4 

242 

19 

37 

37.021 

166.57500 

-0.78100000 

0.00000000 

0 

5 

242 

19 

38 

37.022 

166.81167 

-0.78333333 

0.00000000 

0 

6 

242 

19 

47 

13.026 

168.98083 

-0.80833333 

0.00000000 

2 

7 

242 

19 

57 

7.031 

171.47250 

-0.84666667 

0.00000000 

1 

8 

242 

20 

3 

7.035 

172.99833 

-0.87500000 

0.00000000 

2 

9 

242 

20 

17 

41.043 

176.66667 

-0.95433333 

0.00000000 

2 

10 

242 

20 

34 

51.007 

180.98750 

-1.0716667 

0.00000000 

1 

11 

242 

20 

42 

36.012 

182.93750 

-1.1326667 

0.00000000 

? 

REPLACE  4 

»>  OBS:  SAT#, DAY, HR, MIN, SEC, RA, DEC, RANGE 
? 

83564  242  19  37  37.0  166.576  -.78  0 


USE 

# 

DAY 

HR 

MIN 

SEC 

RA (DEG) 

DEC (DEG) 

RANGE (KM  OR  ER) 

2 

1 

242 

19 

34 

43.020 

165.84583 

-0.77400000 

0.00000000 

1 

2 

242 

19 

35 

43.020 

166  .09583 

-0.77666667 

0.00000000 

1 

3 

242 

19 

36 

37.021 

166.32500 

-0.77833333 

0.00000000 

0 

4 

242 

19 

37 

36.999 

166.57600 

-0.78000000 

0.00000000 

0 

5 

242 

19 

38 

37.022 

166.81167 

-0.78333333 

0.00000000 

0 

6 

242 

19 

47 

13.026 

168.98083 

-0.80833333 

0.00000000 

2 

7 

242 

19 

57 

7.031 

171.47250 

-0.84666667 

0.00000000 

1 

8 

242 

20 

3 

7.035 

172.99833 

-0.87500000 

0.00000000 

2 

9 

242 

20 

17 

41.043 

176.66667 

-0.9 5433333 

0.00000000 

2 

10 

242 

20 

34 

51.007 

180.98750 

-1.0716667 

0.00000000 

1 

11 

242 

20 

42 

36.012 

182.93570 

-1.1326667 

0.00000000 

? 

FILE 

SECURITY  CLASSIFICATION  ? 
U 
? 

EXIT 


Fig. 3.  IOD  — console  session,  part  3. 
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I 

either  of  the  orbit  estimation  methods  via  the  XEQ  command. 
This  is  of  use  in  order  to  check  the  results  of  one  technique 
by  direct  comparison  to  the  other  using  the  same  input  data. 
Any  element  set  estimate  can  be  filed  on  the  Master  Element 
File  (MEF)  via  the  FILe  command. 


III.  COMMAND  STRUCTURE 
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The  program  is  command  structured  with  free  format  input. 
This  means  that  you  tell  the  program  the  name  of  the  function 
you  wish  it  to  perform  and  enough  data  to  perform  that  function. 
Typing  may  begin  anywhere  on  the  line.  Numeric  data  need  not 
be  in  specific  fields.  All  commands  may  be  abbreviated  by  their 
first  three  characters.  Characters  after  the  third  are  ignored. 
Sixteen  commands  are  available.  If  an  unrecognized  command  is 
input  the  program  responds  by  listing  the  set  of  available  com- 
mand abbreviations  on  the  CRT.  Table  1 gives  a graphic  presenta- 
tion of  the  available  commands  and  their  relationship  to  each 
other. 

The  form  of  a typical  command  line  is: 

COM  argl , arg2 , . . . , argn 

As  noted  above,  the  command  field,  "COM"  in  the  example,  must 
be  at  least  three  characters  long.  The  number  of  argument  fields 
depends  upon  the  particular  command  being  issued.  The  command 
and  the  arguments  may  be  separated  by  commas,  equal  signs, 
slashes  or  one  or  more  spaces.  A space  cannot  be  used  in  combi- 
nation with  one  of  the  other  seperators  because  the  operating 
system  interprets  this  as  a null  (zero)  argument. 

A typical  set  of  input  commands  might  look  like  this: 

OBS  83598  242  03 

SELECT 

FIL  95001 


TABLE 


FILe  <1t>  file  orbital  elements  on 


This  retrieves  up  to  12  hours  of  system  observation  file  data 
(or  20  observations,  whichever  comes  first)  beginning  at  day  242 
0300  GMT,  selects  a consistent  and  time  optomized  set  of  obser- 
vations from  the  buffer  and  executes  the  nominal  estimation 
method,  then  files  the  set  on  the  MEF  as  95001.  This  is  a very 
basic  operation  session.  Many  more  features  are  available,  as 
described  below. 


IV.  COMMANDS 

Required  inputs  are  noted  in  upper  case.  Optional  input 
and  argument  paraforms  are  in  lower  case.  Optional  arguments  are 
enclosed  in  left  and  right  carets. 

IV. 1 Satellite  Number 
SATellite  <number> 

If  no  argument  is  input,  the  current  satellite  number  is  listed. 
Number  should  be  between  1 and  99999.  Note  that  the  number  is 
set  to  that  of  the  observations  read  each  time  the  OBS  command 
is  issued. 

IV. 2 Filing  Elements  on  the  MEF 


FILe  <sat#> 

If  no  argument  is  input,  the  element  set  is  filed  under  the 
current  number.  Otherwise,  it  is  filed  under  the  input  number. 
Existing  MEF  records  are  not  automatically  overwritten.  The 
user  is  advised  that  the  condition  exists  and  given  the  option 
of  overwriting. 

IV. 3 Sensor  Coordinate  Specification 
SITe  <name<lat  long  alt>> 

If  no  arguments  are  input,  the  current  sensor  coordinates  are 
listed.  The  user  may  select  any  of  the  ten  internally  stored 
sensor  locations  by  name  (see  below) . If  an  unknown  site  name 
(and  no  coordinates)  is  input,  the  user  is  given  the  opportunity 
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to  input  the  corresponding  site  location  in  double  precision. 

Or,  he  may  input  the  location  on  the  command  line  in  single 
precision. 

The  input  arguments  are: 
name  - 4 character  ASCii 
lat  - degrees  north  geodetic  latitude 
long  - degrees  east  longitude 
alt  - meters  height  above  the  mean  geoid 
The  available  site  list  is: 

MIL  - Millstone  Hill  L-band  radar,  Westford,  MA 

HAY  - Haystack  Observatory  X-band  radar,  Westford,  MA 

ARE  - National  Astronomical  Observatory  UHF  Radar,  Arecibo,  PR 

EGL  - FPS-85  phased  array  radar,  Eglin  AFB,  FL 

STAL  - Lincoln  ETS  optical  site.  Stallion  Range  Center,  WSMR,  NM 

SHEM  - Shemya  AFS  radar,  Shemya,  Alaska 

AMOS  - Optical  site,  Maui,  HI 

EDW  - Baker-Nunn  camera,  Edwards  AFB,  CA 

NAV  - NAVSPASUR  USN  space  surveillence  fence,  Dahlgren,  VA 
RML  - Range  Measurements  Laboratory  optical  site  near  Patrick 
AFB,  FL 

If  a name  of  LIST  is  input  the  list  of  site  abbreviations  is 
written  on  the  CRT. 

IV.  4 Reading  Observational  Data 

The  command  OBS  is  used  to  read  observational  data  into  the 


program  in  one  of  two  formats.  The  choice  of  format  is  deter- 
mined by  the  input  arguments. 


Form  1 reads  data  from  user-prepared  data  (from  the  card 
reader  or  terminal  keyboard) . 

OBS  unit  <datatype> 

unit=ADDS  , read  from  the  ADDS  CRT  keyboard 
unit=CR,  read  from  the  card  reader 

datatype=RA,  read  right  ascension/declination  data 
datatype=any  other  value,  read  azimuth/elevation  data. 

The  order  of  the  input  is  satellite  number,  day,  hour, 
minute,  second  (GMT) , RA  (or  AZ)  DEC  (or  EL) , range.  Angles 
are  expected  in  degrees.  Range  may  be  a KM  or  earth  radii, 
as  you  choose.  Data  is  read  in  free  format  and  double  preci- 
sion, prompted  at  the  terminal.  Range  should  be  supplied  as 
zero  for  angles-only  data.  The  observation  stream  is  terminated 
by  a $$  or  by  a change  in  satellite  number.  Note  that  this 
means  one  observation  is  lost  if  observations  on  different 
satellites  are  juxtaposed  without  an  intervening  $$. 

Form  2 reads  recorded  positional  data  from  the  system  obser- 
vation file. 

OBS  sat#<day<hr<min<sec>>>> 
sat#  is  the  desired  satellite  number 

day  hr  min  sec  is  the  GMT  time  (as  accurately  as  you  need  to 
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enter  it)  of  the  first  observation.  Up  to  20  observation  sets 
or  12  hours  of  data  (whichever  occurs  first)  will  be  retrieved. 
The  time  of  each  observation  found  subsequent  to  the  input  start 
time  will  be  listed  on  the  CRT  screen.  For  either  form,  the 
observation  buffer  is  time  sorted  and  listed  on  the  CRT  after 
the  last  observation  is  read. 

IV. 5 Adding  Observations 
ADD  unit 

unit  is  described  above. 

This  command  re-enters  Form  1 of  the  OBS  command,  above. 

The  data  type  is  taken  to  be  the  same  as  that  of  the  observations 
already  in  the  buffer. 

IV. 6 Deleting  Observations 
DELete  linel  <line2> 

linel  is  the  first  observation  line  to  be  deleted 

line2  is  the  optional  last  line  to  be  deleted 

The  buffer  is  compressed  and  listed  after  a delete  request. 

IV. 7 Replacing  an  Observation 

REPlace  line 

line  is  the  line  number  to  be  replaced 

A separate  read  in  Form  1 is  queued  to  the  terminal  to 
actually  replace  the  data.  The  buffer  is  time  sorted  and  listed 
after  each  REPlace  request.  The  use  code  associated  with  each 
observation,  indicating  the  number  of  times  each  observation  is 
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the  object  of  selection,  is  set  to  zero  for  the  line  replaced. 

IV. 8 Editing  Angles  and  Range 

R. A . line<angle> 

DEC  line<angle> 

RANge  line<distance> 

line  is  the  line  number  of  the  observation  to  be  edited. 

If  no  other  argument  is  input,  a seperate  double  precision  read 
is  queued. 

angle  is  expected  in  degrees. 

distance  is  expected  in  either  KM  or  earth  radii. 

The  use  code  for  edited  lines  is  set  to  zero.  The  buffer 
is  listed  after  each  edit  reqeuest. 

IV. 9 Listing  the  Observations 

List  <linel<line2>> 

If  no  arguments  are  input,  all  observations  are  listed. 
linel  is  the  optional  starting  line  number 
line2  is  the  optional  ending  line  number 
The  listing  is  written  on  the  terminal  only. 

IV. 10  Observation  Selection  and  Orbit  Estimation 

SELect  <i,j<k>> 

If  no  arguments  are  input  an  automatic  selection  method  is 
used.  i,j,k  are  the  optional  line  numbers  of  the  observations 
to  be  chosen.  Three  observations  are  necessary  for  range  esti- 
mation and  for  the  Herrick-Gibbs  orbit  estimation  method.  Two 


angles  plus  range  sets  are  needed  for  the  true  anomaly  method. 
The  use  code  is  incremented  each  time  an  observation  is  selected 
Orbit  estimation  is  performed  after  the  SELect  command  is  issued 
IV. 11  Orbit  Estimation  From  Selected  Observations 
XEQ  <method> 

This  command  is  issued  to  generate  elements  via  the  alter- 
nate estimation  method,  or  re-generate  the  elements  from  the 
primary  method  without  reselecting. 

method=HG  to  execute  the  Herrick-Gibbs  method,  otherwise 
execute  the  true  anomaly  method. 

The  case  number  is  not  incremented,  nor  are  the  use  codes. 

IV. 12  Minimum  Estimated  Error 

MEError 

This  command  lists  the  case  number  having  the  smallest 
true  anomaly  iteration  convergence  parameter  (and  places  that 
element  set  in  a position  from  which  it  can  be  filed  on  the  MEF) 
Though  this  does  not  guarantee  the  quality  of  the  orbital  ele- 
ment set  for  that  case,  it  is  an  indicator  that  these  orbital 
elements  (and  those  with  similar  convergence  parameter  values) 
are  reasonable.  The  CASE  numbering  and  Minimum  Error  Estimation 
are  begun  anew  each  time  the  OBS  command  is  issued.  The  MEE 
command  lists  its  results  on  both  the  CRT  and  printer  file. 

IV. 13  Specifying  Epoch  Year 


YEAr  < number > 


V.  USER'S  NOTES 

V.l  Structure 

Certain  commands  are  linked  together,  so  that  one  is  always 
issued  by  the  program  after  the  user  has  issued  another.  In 
particular,  the  SITe  command  is  issued  by  the  program  when  it 
begins  execution  in  order  to  set  the  default  site  coordinates 
(STAL) . List  is  issued  by  the  program  after  any  observation 
editing  (OBS/ADD/DEL/REP/R. A. /DEC/RAN) . XEQ  is  issued  by  the 
program  after  the  SELect  command  is  issued  by  the  user. 

Both  SEL  and  XEQ  check  to  see  that  sufficient  data  is  avail- 
able to  execute  the  desired  estimation  routine.  If  not,  a mes- 
sage is  written  on  the  terminal  and  the  next  command  is  prompted. 

V.2  Limitations 

Because  the  program  has  little  information  about  the  nature 
of  satellite (s]  whose  orbit (s)  is (are)  to  be  estimated  it  can 
do  little  to  insure  the  integrity  of  the  observations.  In 
particular,  the  user  is  cautioned  against  loading  data  which 
spans  more  than  one  satellite  orbit  (or  pass,  whichever  is  less); 
from  more  than  one  sensor;  or  on  more  than  one  satellite  in  the 
same  observation  read  request.  The  result  of  violating  these 
admonitions  could  be  disastrous. 

V. 3 Automatic  Selection  Algorithm 

The  automatic  selection  algorithm  uses  the  estimated  cir- 
cular mean  motion  for  RA/DEC  data,  or  1 rev/day  for  AZ/EL  data. 
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It  tries  to  locate  data  spanning  1 hour  divided  by  the  mean 
motion  estimate.  The  third  observation  is  not  allowed  to  be 
closer  than  half  this  span  to  the  first;  the  second  must  be  at 
least  .1  of  this  span  from  the  first.  At  the  same  time,  obser- 
vations with  a low  use  code  are  preferred.  This  prevents  the 
same  sets  of  data  from  being  selected  each  time.  For  ten  obser- 
vations, evenly  spaced  over  one  hour  for  a synchronous  satel- 
lite, 15-20  different  observation  selections  can  be  made  before 
the  process  begins  to  repeat  selections.  By  that  point  it  should 
be  obvious  that  either  a good  orbit  has  been  generated  or  more 
data  is  necessary. 

V. 4 Selecting  a High  Quality  Orbital  Element  Set 

No  estimation  technique  will  produce  high  quality  values 
for  the  mean  motion.  It  is  a well  known  fact  of  sampling  theory 
that  if  one  wishes  to  estimate  the  period  (mean  motion)  of  a 
periodic  function  (orbit)  he  must  observe  it  over  several  repe- 
titions (revolutions) . So,  with  only  5-10%  of  an  orbit  we  can 
make  only  a rough  estimate  of  its  mean  motion.  Better  estima- 
tion could  be  done  through  a differential  correction  process. 

For  orbits  of  eccentricity  .1  or  less,  the  estimated  circular 
mean  motion  (which  is  really  the  apparent  angular  velocity  over 
the  observation  span)  is  probably  a better  value  for  mean  motion 
than  that  computed  by  the  orbit  estimator.  This  value  can  be 
substituted  for  the  computed  value  using  the  satellite  batch 
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program  SAT  (or  RSAT  for  the  real-time  operator) . 

Judging  the  quality  of  orbit  estimates  is  not  easy.  The 

MEE  command  provides  some  measure  of  quality  assurance.  The 
best  indicator  of  quality,  however,  is  consistency.  Check  the 
convergence  parameters  to  see  that  they  are  not  large.  Check 
the  inclination,  right  ascension  of  the  ascending  node  and 
eccentricity  for  consistency.  The  other  two  angles  (mean 
anomaly  and  argument  of  perigee)  can  change  dramatically  for 
low  inclination,  low  eccentricity  satellites.  For  higher 
inclination  (above  10°),  higher  eccentricity  (above  .2)  satel- 
lites, the  argument  of  perigee  can  also  be  checked  for  consis- 
te  ncy . 

Grossly  different  values  of  mean  motion  tend  to  indicate 
poor  elements.  Table  2 summarizes  these  checks.  In  general, 
the  longer  the  observation  span  the  better  the  mean  motion 
estimate  and  the  poorer  the  plane  (inclination,  R. A.  of  node 
and  argument  of  perigee) . The  latter  is  because  after  some 
time,  perturbations  to  the  orbit,  which  are  ignored  in  the 
estimate,  become  significant. 
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VI.  PROGRAMMER'S  NOTES 

IOMOD  is  an  overlay  structured  module  generated  by  M4EDIT . 

It  contains  nine  level  1 segments  and  the  main  program.  The 
level  1 segments  are  subroutines: 

ALTRES  alternate  resource  assignments 

IOOBS  observation  input  and  editing 

OBSEL  observation  selection 

SITPAR  site  location  selection 

RANGE  range  estimation 

IONU  true  anomaly  orbit  estimation 

I0HG3  Herrick-Gibbs  3 observation  orbit  estimation 

PFUEL  element  set  conversion  to  classical  and  printing 

IOFILE  element  set  fileup  on  the  MEF 

The  module  is  cataloged  on  disk  partition  TM  as  IOD.  RIOD 
is  a simple  duplicate  to  IOD.  Subroutine  ALTRES  is  called  when 
program  execution  begins.  It  is  told,  "If  I was  activated  as 
RIOD,  assign  my  terminal  resources  to  the  ADDS  symbiont."  This 
is  a slight  over  simplification,  but  conveys  the  intent  of  the 
call.  In  testing  it  was  found  that  the  simple  duplicate  did  not 
follow  the  module  if  it  were  recataloged.  Hence,  if  recataloging 
is  necessary  the  programmer  is  advised  to  reduplicate  as  well. 
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APPENDIX 


SAMPLE  RESULTS 

Sample  estimates  were  produced  for  four  satellites  in  three 
orbit  classes.  For  each  satellite  the  estimates  were  made  with 
and  without  range  for  the  same  set  of  angle  data. 

The  data  used  were  generated  by  propagating  the  "actual" 
orbital  elements  and  producing  the  look  points  for  the  ETS  over 
a four  hour  time  interval.  The  look  points  were  printed  to  a 
precision  of  .01  degree  in  angle  and  .1  kilometer  in  range. 

This  is  then  the  accuracy  of  the  input  data. 

The  satellites  used  are  8195  and  8187,  both  high  inclina- 
tion, high  eccentricity,  half  day  period  satellites;  83561,  a 
low  inclination,  low  eccentricity,  one-day  period  satellite; 
and  8820,  a high  inclination,  low  eccentricity,  225  minute  period 
satellite.  Data  on  8195  is  as  it  approached  perigee.  Data  on 
8187  is  as  it  approached  apogee. 

For  most  sample  results  it  is  easy  to  see  that  having  range 
data  makes  a substantial  improvement  in  the  estimated  orbit. 

The  curiously  "good"  results  for  8820  and  the  "inconsistent" 
results  for  8195  angles-only  cases  can  be  shown  to  be  within  the 
error  bounds  dictated  by  the  accuracy  of  the  input  data.  Range 
estimates  are  off  by  2%-3%  for  these  cases. 

The  cases  labled  "A"  were  generated  by  the  true  anomaly 
iteration  method;  those  labled  "B"  were  generated  by  the  Herrick- 
Gibbs  method. 
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